Correlations in a Generalized Elastic Model: Fractional Langevin Equation Approach 
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The Generahzed Elastic Model (GEM) provides the evolution equation which governs the stochas- 
tic motion of several many-body systems in nature, such as polymers, membranes, growing interfaces. 
On the other hand a probe [tracer) particle in these systems performs a fractional Brownian motion 
due to the spatial interactions with the other system's components. The tracer's anomalous dynam- 
ics can be described by a Fractional Langevin Equation (FLE) with a space-time correlated noise. 
We demonstrate that the description given in terms of GEM coincides with that furnished by the 
relative FLE, by showing that the correlation functions of the stochastic field obtained within the 
FLE framework agree to the corresponding quantities calculated from the GEM. Furthermore we 
show that the Fox //-function formalism appears to be very convenient to describe the correlation 
properties within the FLE approach. 



I. INTRODUCTION 

Polymers [J-yl, elastic chains and membranes [3, |j-0] , 
rough surfaces |8l-[l0| can be described by a continuum 
elastic model which accounts for their general stochastic 
behavior. This model, named Generalized Elastic Model 
(GEM), is defined by the following stochastic linear dif- 



ferential equation in partial derivatives 11 1 



j^W,t)+ri{x,t). 



|-h(f,t)= / d^x'Kix-x') 

(1) 

The stochastic field h is D-dimensional and defined in the 
d— dimensional infinite space. The white noise appearing 
in ^ satisfies the fluctuation-dissipation (FD) relation, 
i.e. 



(77j (a-, t) 7?fe (f , t')) = 2kBTA (f ~ x') 6, k6{t - t'). (2) 

{j,k £ [1,-D]) where A (r) corresponds to the hydrody- 
namic friction kernel which couples different sites in x 
and x' through a fluid-mediated interaction. The inter- 
nal elastic coupling is instead embodied by the fractional 
derivative defined via its Fourier transform by [l2| 



•^^^ 19 If I 



^-w 



(3) 



Another common definition is given in term of the Lapla- 



cian A as hi 



d\x\ 



-=-(-A) 



■./2 



The systems whose the dynamical behavior is de- 
scribed by the GEM ([Ij can be divided into two different 
classes, according to the type of hydrodynamic interac- 
tions that characterize them: long ranged or local. 



A. Long ranged hydrodynamic interactions 

This situation occurs when the friction kernel is defined 
by a general expression like 



K{r) 



H 



(4) 



where a <C |r| <C i, i corresponds, for example, to the 
screening length or the maximum size of the system, and 
a is the smallest length scale up to which the continuum 
description furnished by ((T]) keeps its validity. Whenever 
our analysis will require a regularization at small and/or 
long distances, the largest and smallest length scales will 
then come into play as the integral's upper and/or lower 
cut-offs. 

The hydrodynamic interactions are often represented by 
the equilibrium average of the Oseen tensor, which in an 
embedding dg-dimensional space (de > 3) reads, accord- 
ing to (HI): A (f) ~ \rf~'^' [ll,[i3. The following systems 
belong to this class. 

- Fluid membranes. The height of a fiuctuating mem- 
brane is represented by the scalar quantity h{x,t), i.e. 
D — I [2, 5-7]. The point x on the planar surface im- 
plies d — 2, which in turn gives d^ ^ D + d = 3 and 

a — 1. Since for small deformations ( V/i ^1) the 
bending free energy is ex (Ah) flSl], z = 4 in ([IJ. 

- Semiflexihle polymers, h stands for the 3 spatial co- 
ordinates of a polymeric segment {head) while x is the 
strand's 1-dimensional internal coordinate [curvilinear 
abscissa): D — 3, d = 1. The embedding dimension 
in this case coincides with D [d^ — 3), yielding a = 1 
0, Q • The bending elastic energ y a ssociated with the 
chain's deformation implies z = 4 [16| as in the previous 
example. 

- Flexible polymers. In this systems h and x still cor- 
respond to the position and the curvilinear ascissa of the 
monomer in the polymeric chain: D — 3, d — 1. However 
the beads interaction is just represented by an harmonic 
coupling [z = 2) and the Zimm's equilibrium approxi- 
mation of the Oseen tensor in 9 solvent gives a — 1/2 



B. Local hydrodynamic interactions 

This kind of systems present no fluid-mediated inter- 
actions, namely A{'r) — S'^{r). This can be attributed 
either to the large screening among the elementary com- 
ponents of the system or to an interaction which is purely 
mechanical. Examples are: 

- Rouse polymers. Here, as in the case of (semi)flexible 
chains, h stands for the bead's 3-dimensional position, 
and X for the bead's position along the chain [l|, [l8| . 
The hookean interaction gives z — 2. 

- Single file system. The system of N hard-core rods 
diffusing on a line without overlapping is known in lit- 
erature as single file (see [I3 and references therein). 
Recently it has been shown that the system's dynamics 
can be reduced, within a very good approximation, to a 
1-dimensional harmonic chain problem {harm,onization) 
where h{x, t) is the position of the x-th particle on the 
substrate at time t. 

- Fluctuating interfaces. In this systems h plays the 
role of a scalar field (mostly the height of a rough sur- 
face in d dimension) which is subjected to a non-standard 
elastic force embodied by the fractional derivative of or- 
der z [1, Q . This is actually the generalization of the 
Edwards- Wilkinson equation for the fluctuating profile 
of a granular surface, for which d = 2, z ^ 2 [4|. In 
systems such as crack propagation [20[ and contact line 
of a liquid meniscus [2l| d = 1 , and the restoring forces 
are characterized by z = 1. 

- Solid surfaces. If /i is a step, namely a line bound- 
ary at which the surface changes height by one or more 
atomic units 10], the value of z in eq.(IT]) is found to 
he z = 2,3 or 4 {d — 1) according to the character of 
the atomic diffusion: periphery, terrace or attachment- 
detachment diffusion respectively. 

- Diffusion-noise equation. In this case h represents 
the density field on a d-dimensional surface x while the 
diffusion operator sets z = 2 [2a |. 

The values of the parameters related to each of the 
models formerly listed are summarized in table ID 

In [ll| we addressed the question of the motion of 
a tracer (probe) particle in the systems whose dynamics 
obeys to ([T|). Although the whole system dynamics is 
Markovian, the particle placed at a position x undergoes 
a subdiffusive motion on the score of persistent memory 
effects due to the spatial correlations with the rest of 
the system. Roughly speaking, a tracer particle experi- 
ences two kinds of interactions: the first one is the cou- 
pling with the surrounding heat bath, whose mathemati- 
cal expression is furnished by the Langevin random force 
ri{x,t) and the corresponding FD relation ^. The sec- 
ond interaction is inherent to the system: the probe par- 
ticle is coupled with the rest of the system through both 
the hydrodynamic term ^ and the fractional Laplacian 
([3]). This "internal" coupling originates the spatial cor- 
relations responsible for the tracer's long-ranged non- 



System 


D 


d 


z 


a 


Fluid membranes [2, ^-_7] 


1 


2 


4 


1 


Semiflexible polymers [2, ^ 


3 


1 


4 


1 


Flexible polymers [1, 17} 


3 


1 


2 


1/2 


Crack propagation [20] 


1 


1 


1 


- 


Liquid meniscus [21] 


1 


1 


1 


- 


Rouse polymers [1, J^ 


3 


1 


2 


- 


Single file systems [19] 


1 


1 


2 


- 


Fluctuating interfaces [4, 8^, 9\ 


1 


any 


any 


- 


Solid surfaces [10| 


1 


1 


2,3,4 


- 


Diffusion-noise equation [22] 


1 


any 


2 


- 



TABLE I. Values of the parameters D, d, z and a charachter- 
izing the GEM ^ for the systems listed in the Introduction. 
Note that for systems whose the hydrodynamic interaction is 
local the value of a is absent. 



Markovian memory effects. On the other hand, the non- 
Markovianity of the probe particle's anomalous dynamics 
leads to the description in terms of fractional Brownian 
motion (FBM) 1231 which obeys a fractional Langevin 
equation (FLE) [l^ |2J|- Within this framework the 
strong internal interactions are mimicked by the colored 
noise term and the fractional derivative, connected by 
the generalized fluctuation-dissipation relation. In this 
article we will show how the representation of the tracer's 
stochastic evolution given in terms of FLE offers the same 
level of accuracy furnished by ([T]). Indeed any kind of 
physical statistical observable can be evaluated starting 
from ^ as well as from the corresponding FLE. 

This paper is outlined as follows: in SecUl] we start 
from the GEM ^ deriving the expressions for the h- 
autocorrelation function (and the corresponding mean 
square displacement) in the case oi z > d, z < d and 
z = d {{d - l)/2 < a < d). In Section [ml start- 
ing from the expression ([T]), we draw the FLE equation 
for the tracer particle placed at position x when z > d. 
In Section IIVI we derive the properties of the noise ap- 
pearing in the FLE and we demonstrate the validity of 
the fluctuation-dissipation relation for the probe particle. 
Section |V] will be devoted to the /i-correlation function 
arised from the FLE framework and we furnish its gen- 
eral expression in terms of the Fox functions. In partic- 
ular we apply the developed formalism to systems such 
as fluid membranes, proteins and fluctuating interfaces, 
recovering results previously derived in literature. In Ap- 
pendix|X]we report the calculations for the hydrodynam- 
ics term in the limiting case a = d. In Appendix [B] we 
provide a demonstration of the appeareance of the Fox 
functions in our analysis, while in Appendix [Cj arc listed 
the main properties of the Fox functions that we make use 
of through our calculations. Lastly, Appendix lEl concerns 
with the derivation of the Generalized Langevin Equa- 
tion for the inter-monomeric distance in a 3-dimensional 
Rouse chain, according to the procedure outlined in [l!|. 



II. ^-AUTOCORRELATION FUNCTION FOR 
THE GENERALIZED ELASTIC PROCESS 

We start from the equation ([T]) . Defining Fourier trans- 
form of tlie stochastic process in space and time as 



{[hj {x, t) — hj {x, 0)] [hj {x, t') — hj (x, 0)]). Because 

of the isotropy of the system henceforth we will drop the 
index j. Since the Fourier transform of the noise corre- 
lation function ^ gets the form {rjj {q,uj)r]k [if ^oj')) = 
2kBT{2Tr)'^+'^A{q}6jkS{uj + uj')S'^{q + q'), after a bit of 
algebra we derive 



h (q, w) 



+00 /'+00 

d 



(fx / dt h (f , t) e 



-i{q-x—ujt) 



(5) 



we find that the general solution of ([ij can be expressed 
in the Fourier space in the following form 



{Sh{x,t)Sh{x,t')) ^ksTx 

-Hoo d''q 9M-t)+g{\i]\.t')-g{\q\.\t-t'\) ,g. 



^ J-00 (2ir)'i 



'Mq)\<l\'t 



h {q, uj) = 



viq,^) 



-iuj + A{q)\q]' 



(6) 



where A (q) is the Fourier transform of hydrodynamic 
friction kernel, 



A(q) 



(47r)''/^ r((d-a)/2) i^g-d ^ y^ i^g-d rf-1 



2-'''"^ln 



r(d/2) 



r(a/2) 



2 

a — d. 



< a < 



When {d — l)/2 < a < d the general expression of (|9]) 
has the following form: 

{Sh {x, t) 5h {x, t')) = K [fit) + f (f) -f{\t- t'\)] . 

(10) 
Three different situations can occur. 

i) z < d {a > d — z). The integrals in © are diver- 
j gents when |(f| — > 00: hence once again we are compelled 
to introduce the cut-off a. Performing the integrations 
yields 



To get eq.© we use the expression for the d-dimensional 
Fourier transform of the isotropic function 0(|r|) [25| 



C:: d'^re-^^'MM) = 

where Jd/2-1 represents the Bessel function of fractional 
order d/2 — 1. The limiting case a = d requires a reg- 
ularization of the integral in the g-space: this has been 
done introducing the cut-off a. Note that if we resort to 
the different definition A (r) ^ ^ ,"'", ^,j , we get the same 

asymptotic expansion (for small \q\) obtained in ([7|) for 
a ~ d\ va this case however, the complete results for 
d = 1, 2 and 3 are presented in Appendix 1X1 The case 
a — {d— l)/2 which requires the infrared cut-off L won't 
be treated in the following. 

In the local hydrodynamic situation A (r) — 5 [r), which 
corresponds to take A (g) = 1 in (jB]) . Therefore the cal- 
culations for local hydrodynamic systems can be either 
performed starting from ([5]) where A (g) — 1, or sim- 
ply setting A = const and a = d va. the correspond- 
ing long-ranged hydrodynamic expressions: this substi- 
tution, which is not to be intended as a limit, allows 
to easily shift from power-law to local hydrodynamics 
throughout the following analysis. Indeed it corresponds 
to a formal procedure to obtain the results for local hy- 
drodynamic systems, starting from the equivalent quan- 
tities elaborated for long ranged hydrodynamic models. 
As a consequence, the case a = d and the ensuing loga- 
rithmic behavior (O, does not have to be confused with 
the systems with local hydrodynamic interactions. 
We then define the ^.-autocorrelation 

function as {dh{x,t) Sh{x,t')) = 



K = 

fit)- 

where 



(2ir)'T(d/2) a+z~d 

a+z — d 



A U^\ 



(-A^(f) 



ct-'rz — d 



/3 = 



t]t^ 
(11) 



(12) 



and the function 7(0, x) is defined in [31.]. It is straight- 
forward to verify that in this case the mean square 
displacement gets to a constant value, {6'^h{x,t)) — >■ 

(2xmJ72)(d-z) (f )''"■ Physically this means that the 
system is overconnected and asymptotically any probe 
remains trapped around its initial position. 

ii) z — d. The integrals in Q exhibit a logarithmic di- 
vergence. After the regularization through the insertion 
of a the result is achieved substituting in PH)) 



2kBT-!T' 



d/2 



K 

(27r)'T(d/2)Q 

/(t) = i?i(A(f)"i)-Hln[A(f)"t]+C, 



(13) 



where C is a constant value that we set to and Ei{x) 
denotes the exponential integral [31]. In this case 
the system attains an asymptotic logarithmic diffusion: 

{5MS^.t)) ^ (4«rwfcln[^(S)"*]- This is a bor- 
derline case, where the probe is not completely free to 
diffuse away from its initial position: the ensuing erratic 
motion is then logarithmic. 

Hi) z > d. This is the most interesting situation: in 
this case the integrals can be performed effortlessly to 
give 



7^ _ 2kBTTT'^''^ A^T{1-P) 
" (27r)<ir(d/2) z-d 



fit) = t^. 



(14) 



Therefore, the tracer particle placed at a given x per- 
forms a subdiffusive fractional Brownian motion (FBM) 
[23] given by {d^h (x, t)) = 2Kt'^ . In the language of fluc- 
tuating interfaces an interface is called rough in this case 

In the following Sections we will focus on the situation 
in) for which z > d. 



III. FRACTIONAL LANGEVIN EQUATION 

Starting from the solution (|6]) we now discuss the 
derivation of the fractional Langevin equation for the 
probe particle placed at a given position x. 
Let us first multiply both sides of ^ by K^{—iuj)^. As 
one can see from the derivation presented below, this 
is the only choice leading to the FLE which obeys the 
FD relation. According to [ll|, choosing another ar- 
bitrary power instead of {—iuj)" would lead to another 
equation among the family of Generalized Fractional 
Langevin Equations (GFLE), but without FD relation 
fulfilcd. Thus eq.([n|) becomes 



K+{-iuj)^h{q,uj) 



-lU} + A\q\ 



K^ is a constant that we introduce in order to fulfill 
the fluctuation-dissipation relation: indeed this physical 
constraint must be satisfied regardless of the description 
made of the tracer's dynamics, i.e. both by the Marko- 
vian Langevin description given in ([1]) and by the frac- 
tional Langevin representation that we are aimed at de- 
riving. 
Eauation (|15p can be rewritten as 



K 



^{~iLo)^h {q,u) = K+r^ {q,Lu) <D (<f,c^) , 



(16) 



where we introduce the function $ (x, t) whose Fourier 
transform in space and time reads 



To transform the left hand side of Ea. p6|) . we introduce 
the Caputo derivative with lower bound equal to — oo, 
which for a "sufficiently well-behaved" function (hit) is 
defined as follows 1131 w3 , 



oUit) = 



1 



dt'- 



r (1 - /3) J_^ (i _ t')^ dt 

and whose Fourier transform reads as 



Vl^'^(''^' o</^<i' 



Aujt T~\P 



dte"^'Dl<j){t) ^ {-iLj)^<l){uj). 



(19) 



(20) 



Thus, we get finally the fractional Laiigevin equation for 
the stochastic field h (x, t) [gI. [TTI. [HIm [26} , 



'+r)/5 



K+D'lh{x,t) ^C{x,t) 



(21) 



From Eq. ([2T|) and the definition (|T9|). the requirement of 
the validity of the fluctuation-dissipation relation reads 



{ax,t)C{xX))=kBT 



K^ 



T{l-P)\t~t'f' 



(22) 



which relates the autocorrelation function of the noise ( 
standing in the right hand side of eg. ([21] )) to the damping 
kernel (determined by the fractional derivative (J19p ) in 
the left hand side of ea. (pT|) . From this requirement the 
value of K'^ will be set. This will be done in the next 
Section. It is important to remark that the random fleld 
(fT8| is still Gaussian, since it is the linear combination of 
the Gaussian noise ?7(x, t) but is also power-law correlated 
in time according to the FD relation (P^. hence is a 
fractional Gaussian noise (fGn). 



IV. FRACTIONAL GAUSSIAN NOISE 
CORRELATION FUNCTION 

We want to evaluate the two-point two-time correlation 
function of the fGn appearing in (PT|) . From ^ and from 
the definition ([T5| we obtain 



$ (q,uj) 



{-iuj)f^ 



-iuj + A\c[\ 



z-\-ct — d ' 



(17) 



Now, we simply make an inverse Fourier transform of 
Eg. dTGl) in space and time. In the right hand side we 
get a new noise term C [x, t), which is the convolution of 
$ (x, t) with the white Gaussian noise 77 (x, t), i.e. 



(:{x,t)^K 



- T<-+ 



d'^x' 



dt'T]{x-x',t-t')<^{x',t'). 

(18) 



where 



95 



{C{x,t)c(^x',t'^) = 2kBTK+^x 



(23) 



(x — x' ,Uj] = 
J dxidx2A {xi ~ X2) ^ {x — xi,uj) ^ ix' ~ X2, —uj] 



(24) 



To calculate (p4)) we make use of the hydrodynamic 
term's Fourier transform ([7]) and of P7|) : after straight- 
forward manipulations it takes the form 



^ (x — x'jOj] = 



(2iry/2 I 



1 2/3 



x/o d\^ 



k1° 



'Jd/2-l(l'?l|^-a:'|) 



l--d/2 



(25) 



which, after a change of variable, becomes 



If [x — a;', 1^ ) 



^2(z + c,-d) 



(27r)<i/2 



/o°°c?y 



-d/2 



1+y 



2TT+^ 



dy^d/ 



. -2c.-d + 2 
^|2(^ + Q-d) 

l/(^+Q 



2-1 



l-d/2 



al'^l' 



j^l/[z + a,~d) 



(26) 

To proceed further, at this stage we employ the formalism 
of the Fox i7— functions. These functions, introduced by 
Fox in 1961 33|, are special functions of a very general 
nature which allow us to present the results in a uni- 
versal and elegant fashion. For a general theory on the 
if— functions we address the reader to the monograph 
of Mathai and Saxena [SJ] and to ref. [S^l . As interest- 
ing applications of iJ— functions we could mention an ex- 
actly solvable model of linear viscoelastic behavior '36] , 
the 7J— function representation of non-Debye relaxation 
|37l ISSj and of the solution of the space-time fractional 
diffusion equations [3^, EO]- We present the defintion 
and the basic properties of the 7f— functions in Appen- 
dices [B] and O respectively. Moreover in Appendix [B] it 



2{z+a-d)\ 



is shown that the function 1/(1-1-?/' 

the expression (|26p can be cast in term of a Fox function: 



appearmg m 



1 



- -^11 



l + y7 J 
where we introduced the short notation 



y 



0,i 



'O'i 



J = 2(z + a-d) 
After this substitution (pS)) reads 



(27) 



(28) 



(fi (x — x' ,Uj] 



-2e.-d+2 



(2^)d/2^ 




l-d/2 



with k = (|w|/A)2/-r 



(29) 
Using the property (|C4p 



the integral can be evaluated to give 



(f Ix ~ x' ,uj] 



fc° + l-<J/2 



frl2 
-"31 



(2^)d/2^|^| ^ 

1-f,^) (o,i) (1 



l-d/2 



a+2-d 1^ 
2 ' 2; 



(30) 



which, thanks to (IC2[) and (jC3p . gets the form 



iy9 I a; — a; , oj 



./3+- 



-2Q-d + 2 



W 



ff2 1 



2(27r)<'/27 

' 4z-|-2Q:-3d-2 j_ 



l-d/2 



l' d-2 j_^ 
V 4 ' 2/ 




' 2-d j_^ 
. 4 ' 2/ 



(31) 

Applying again the property (|C3p . the noise correlation 
function (1231) can be written as 



2kBTK+^A 



+ 2 .2P-1 



{C{x,t)C(x',t'))- ,. + d^^d/2+l 

xcos(a;(i-i'))ii2i 



a: — x' 



/o ^^'^^ 



( ~a 1\ a I] ( 2-a-d 1\ 
V2'2/ V^'7y \ 2 '2/ 

(32)" 

where /3 and 7 are given by ((T^ and (|28p respectively. 
The integral in (|32|) can be solved by referring to the 
property (ICSp of the Fox function: the final expression 
for the noise correlation function then reads 



(C(x,t)c(x',t')) 



ksTK+^A 



+2 42/3-1 S-o;' 



xHll 



A\t-t'\ 



2c, + djr(d+i)/2 

'1 1 
>2' 2 



■J'ZJ V^' 2 j 



^ 2-et-d 2^ 
2 ' 4J 

(33) 

The former expression is the central result of this paper. 
It states that the fGn entering the FLE (|2ip is not only 
correlated in time, as required by the physical constraint 
(l22t . but it is also space-correlated. This means that the 
space correlations appearing in ([IJ, which are embodied 
by the hydrodynamic term as well as by the fractional 
Laplacian, are translated into space-time correlations of 
the noise in the FLE dynamical representation. 

We recall now that the coefficient K'^ in (I55|) and (P^ 
is still undefined, our aim is to set its expression. For 
this purpose we calculate the autocorrelation function of 
the noise, i.e. we set a? = a:' in p3l) . To do so, we first 
expand the Fox function for small argument according 
to (|C7p and restric ourselves to the main term which 
diverges at x = x': 



tt22 

-^33 



A|t-t'| 



(i^) iP^h) (0.^) 



2^+°V¥A^-'^ r(/j) \t-t'\ 

7 r(d/2) \3-x'\°' ■ 

(34) 
We then plug such expression in ([33]) . achieving the fol- 
lowing final form for the fractional Gaussian noise auto- 
correlation function : 



{C{x,t)Q{x,t')) 



^bT _^+2 



\t-t'\ 



'll^' 



77^/27 r(d/2)' 



(35) 



Thus the comparison between (|35l) and ((22|) yields the 
value of K^ , namely 



7^+ = (47r)'^/2-i7sin(^/3)^^. (36) 



TWO-POINT TWO-TIME ^-CORRELATION 
FUNCTION 



It follows from (gO]) and (|39l) that the integrand in ([39)1 
diverges as ^ li^l" ■ Therefore, in order be measur- 

able, any physical observable has to be arranged in such 
a way that the divergence will be eliminated. In the fol- 
lowing subsections we will analyze some specific case of 
statistical quantities constructed as linear combination of 
(j39l) and compare them with the corresponding quanti- 
ties arising from ([6l). 



In this Section we address the problem of the evalua- 
tion of the ^.-correlation function within the framework 
of the FLE. Indeed, any kind of statistical observable can 
be expressed in term of correlation functions, whose ana- 
lytical expression can be derived either starting from the 
eq. ([1]) or from eq. ([21]) . Although the representation of the 
system dynamics can be different, the correlation func- 
tions must coincide, since the observable physical prop- 
erties do not have to depend on the chosen description. 
Therefore we can furnish the general expression of the 
correlation function starting from the solution of (l?T|) . 
The solution of (|2ip in the Fourier-Fourier space reads 



A. Mean Square Displacement 

The mean square displacement of the probe particle 
has been defined in secton [ll] as the limit i — >■ i' of the 
autocorrelation function Q ; in term of elementary com- 
ponent p9|) it can be expressed as 



{5'^h (x, t)) = 2 [{h^ (f , 0)) - {h [x, t) h (x, 0))] . (41) 

Using the expansion (pHl) the expression (HTI) can be 
casted as 



h (g, w) 



C(g',^) 

A'+(-^CJ)'3' 



(37) {5^h{x,t)) 



and consequently the two-point two-time correlation 
function is 



{h{x,t)h(x',t'\) 



1 



duj duj' 



e-''"*e-*'^ 



,^,(C(f,c.)c(:r'',c.')) (38) 



^+2 7 47r2 ^ ^ {-iujf{-iuj'f 

Using (|32|) . the Fox function expression for (|38l) reads 



{h{x,t)h(x',t'\) 

s. ™K"(*-*')) u2 1 
^ „23 -«13 



2a + d^^d/2 + l 



/o 



+ 00 



duj'x 






2-Q-d 1 



(39) 

The previous expression constitutes the elementary com- 
ponent starting with whom, any kind of physical ob- 
servable is constructed. Nonetheless, the integral ap- 
pearing in p9|) cannot be solved as we performed in 
([5^ . The reason is that such integral is divergent in 
the limit lo — ^ 0, as it is apparent by an expansion of 
the Fox function at small argument, i.e., recalling that 



tt2 1 



2/7 



2ir 



X — X' 



and using (IC7|) . one has 



[P^h. 



{^^\) /^.^ ( 



r(d/2)cos(^/2) 



2~a-d l\ 
2 ' 2) 
1-13 



(40) 



|a.| 



3T{An)-''/^Al^ /•+°° , l-cos(wt) 

CliJ — 



jrid/2)cos{n(3/2)Ja 



H 



1+/3 



(42) 

Solving the integral we get the expression (I14p obtained 
previously from the GEM dH). 



B. Dynamic structure factor of fluid membranes 

In ref. Q Zilman and Granek derived the short length 
and short time behavior of the dynamic structure factor 
of the fluid membranes. As mentioned in the Introduc- 
tion, the fluid membrane dynamics correspond to take 
_D = 1, d = 2 and z = 4 in ([Ij, moreover the hydrody- 
namic friction kernel is expressed as 



A(r) 



1 



Stt^H 



(43) 



which gives a = 1 in (j4]), ^ is the solvent viscosity. Note 
that unlike the original model we set the bending mod- 
ulus K = 1. According to ([7]) and the definition (j43|) the 
constant A is found to be ^ = 1/(4^). 
The quantity which has been studied in [Ts] is the two- 
point correlation function 



{ h{x,t)-h(x',0\ ) = 2 {h'^{x,0)) -{h{x,t)h(x',0\) 

(44) 
which is shown to be linked to the membrane's dynamic 
structure factor. The last expression can thus be recast 
in term of the quantity ([M)) : using the expansion PO)) 
and the numerical values of d, z, a and A we find 



h(x,t)~h(x',0) ) = 

/+°° rfw I ^ - cos (ujt) i^^i^ 



2kBT 
37r(44)2/3 

X HU 



(4C^) 



'2 1^ 
k3' 6^ 
-1 1\ /'2 1^ 
'J V3' 6^ 



\3' 6/ 



1 1^ 
2' 2y 



(45) 

After changing variable {y = x — x' (4^0;)"'^' '^/2) and 

making use of the property (jC3l) the following general 
simpler expression is achieved for the two-point correla- 
tion function: 



{h{x,t)-h{x'M )^^Io dy 



Hi 



(0,^; 



(1 1^ 

V6' 61 
(5 1^ 
V6' 6J 



(0,^; 



(46) 

Although eq. pS)) furnishes a compact analitycal expres- 
sion for the correlation function at any time and any 

distance x — x' , the integral cannot be computed es- 

plicitely since it displays a logarithmic divergence in the 
limit y -^ 0. Again this can be seen by expanding the 
cosine and the Fox function in (1461) to the second order: 



Hf 



21 



'5 
.6' 6) 

V^' 2 J V6' 6^ 



(5 1) 
V6' 6J 



(o^i] 



47r-27ry^+7ry^. (47) 



We first analyze the limit t = which corresponds to 
the static correlator ( h{x,0) — h Ix' ,0j ), describing 

the membrane's roughness. In this case we immediately 
get from pS)) 



h{x,0)^h{xl0)] ).i^/;-d,^ 



1 rr21 

477'"13 



(0,^) (0) (o,i) 



(48) 



On the score of the previous discussion, using the ex- 
pansion (I47p the above expression can be approximated 
as 



( h{x,0)-h(x',0\ 



( r\ 

\J \x — x' 



4:1T 



dy 



iSldyy 



(49) 



where L stands for the long scale cut-off representing the 
membrane's size. The underlying assumption in (j49p is 
that the major contribution to the integral in (|48p comes 
from y < 1; this, in turn, justifies the lower cut-off ap- 
pearing in the first integral of (|49p : indeed the mini- 
mum relaxation frequency of the membrane's bending 



modes will be given by ujq 



to yo 



■J-) jF which corresponds 



/ L . Hence solving P^ we obtain the 

expression which coincides with that found in [7] for the 
static correlator, i.e. 






valid whenever 



( h{x,0)-h(x',()\ ) 
In 



x — X' 



0.25 



(50) 



< L. The value of -0.25 of the 

correction term can be improved in a regular way; see 
Appendix |DJ 

For the dynamic correlator we can still consider the main 
contribution to the integral arising from y < 1 and use 
(H71) . then 



h{x,t)-h[x',Q)\ )~i^ 

l-cosi I -^3L 



/o dy- 



lo dy- 



1 fi 

2 



(51)" 



Changing the variable back to oj, for times in the intcr- 

3 

"^"^ ' ^ < i < 4^ (-1)^ we can safely 



mediate range 4^ 



replace the upper bound of the integrals to oo achieving 



h{x,t)-h[x',0 )^^ 



knT 



(4?) 



-2/3 



(52l 

Thus we can compute the first integral and introduce the 
infrafred cut-off loq to regularize the second, therefore the 
final expression for the correlation function reads 



( h{x,t)-h[x',{)\ ) 
r(i/3) t'/'' '^^(' 



2tt 



(53) 



(4?)=/^ \S-xf 



where Ci is the cosine integral |31|]. The result ^ 
matches and amends the corresponding expression fur 
nished in Ml. 



C. Donor-Acceptor correlation function in proteins 

In Refs [2^ [23] it has been shown that the dynamics 
of the donor-acceptor (D-A) distance within a protein 
can be mapped into the motion of a fictitious particle 
obeying a FLE with fractional derivative of order 1/2, 
in presence of an harmonic potential whose frequency Wg 



could be phenomenologically inferred a posteriori from 
the experimental data. The detected quantity was the 
autocorrelation function of the D-A distance A.o~A{t) 
that was shown to display an asymptotic Mittag-Lefher 
decay in accordance with the FLE prescription. In or- 
der torecover the experimental results, in Refs [28|, 123 
and |3Q] the authors used a respectively continuous and 
discrete Rouse model accounting for the protein confor- 
mational dynamics: this, in turn, corresponds to take 
£) = 3, d = 1, z = 2 and A (x - x') = J'' (x - x') in (P) 
(see Introduction). The simple Rouse model was shown 
to reproduce the Mittag-LefHer decay of the ^d-a au- 
tocorrelation function for large t and, on the other hand, 
it was shown to lead to the correct 1/2-FLE for the D- 
A distance within the framework developed in [ll|, [l^] , 
with the frequency Wq oc 1/ (xa — xd)- We now want to 
calculate the D-A autocorrelation function without re- 
sorting to the derivation of the FLE for the D-A distance 
(which will be done in Appendix [E|) . but starting from 
the FLE (PT|) valid for the monomer placed at position x 
measured along the protein's profile. 

According to [2^ we define the ^d-a autocorrelation 
function as 



^(,) ^ 2^/^7n-(.)[sin(^)+cos(^)] ^^^^ 



which can be inverted [425 to give 



f72 



V4' 4/ 
v'-'' 2) V4' 1) V2' 2) 



(59) 



2^/^^e-^y (sin ^/2y + cos y/2y) . 
Plugging ([5^ in ([55]) and integrating by parts, we obtain 

{^D-A{t) ■ ^D-A{t')) = 2^^^ \XA - Xd\ X 

X [+2C/o°° dysin (Cy^) g-^w (sin\/2y + cos %/2y) + 

+2V2 /„°° dy '^-^^'-^^'^K -^^y - 2C/„°° dysinCy^ ] 

(60) 

where C — -. — ■ U. Let's analyze the three integrals 

\xA-xnr ° 

within the square brackets. First we consider the follow- 
ing general expression 



{/^D-A{t) ■ AD-A{t')) = 3{AD-A{t)AD-A{t')) (54) 

which can be rewritten in term of the spatial positions 
of the donor and acceptor sites, h{x]j,t) and h{xA,t) 
respectively, as 



dyye-^y sin (Cy^) sin {(3y) = ^/^^e^^. (61) 



We then integrate with respect to /3 both the RHS and 
the LHS of (ETl) to achieve [il 



{AD-A{t) ■ AD-A{t')) = 

6 {{h {xA,t) h {xA, t')) - {h (xA, t) h ixD,t'))) 



(55) 



Therefore, putting in ([39]) the numerical value of the pa- 
rameters and changing variable {y = \xa — X]j\ y/\uj\/2) 
we get 



n e 2U 



dye-^y sin {Cy') (sin (/?y) + cos (/3y)) = ^ ^ ^3/2 

(62) 
Multiplying both sides by 2C and setting (3 — \f2 eq. (|62|) 
yields 



(Az5_^(t) • /^D-A(t')) - i^^^ |a;^ _ ^^1 



C+oo J 
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\^A '•'d\ 
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1 jj2 1 
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V*^' 21 V4' 4/ V2' 2) 



(56) 
Expanding the 7J-function in the former expression gives 



tt2\ 
-^13 



(o,i) 
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V4' 4/ 

(3 l^ 

V4' 4/ 



'1 1^ 
^2' 2; 



2?l^^{\-2y^\ 



(57) 

which guarantees that the integral in ([5^ does not 
present any divergence: we then seek for its explicit so- 
lution. First we express the Fox function's Mellin trans- 
form following the definition (jBll) : according to (|B2I) it 
is found to be 



2(7/?° dysm {Cy^) e"^^ (sin (72?/) -f cos (v^y)) 



^e-^. 



(63) 

We next consider the second integral in (|5(I)) . hereby 
named l2{C), and differentiate it with respect to C, i.e. 



^/2(C)--^^^ 



(64) 



The former differential equation can be solved thanks to 
the initial condition 1431 



W) = 2V2rdy'-^^^y^e-^^y = ^ (65) 

Jo y v2 



yielding 



h{C) 



V2 



erf 



(66) 



where erf denotes the error function [3l| . 

The third integral in (l60t can be evaluated effortlessly 



2C 



dy sin Cy^ = — 



(67) 



Putting dMI, dini) and §7^1 in §^ and substituting the 
expression of C, the final form of the D-A autocorrelation 
function is achieved: 



(Ao-Ait) ■ AD-A{t')) 



SVIfcflT 



XA- Xd\ X 



2TT\t~t'\ 



\xa~xo\ 



4|t-t'| 



I TT r \ \XA-Xo\ / 1 11 



which is exactly the expression found in ref [19]. More- 
over (p5|) recovers the asymptotic decay of the autocorre- 
lation function Cx (t) observed in the experiments [2a, UM 
that indeed was found to be 



C'x(t) — —El/2,1 



knT 



OJn 



''-f, (69) 



where wg is the characteristic frequency of the potential 
and to — l^j , with ^ generalized damping coefficient. 

On the other hand the expression (j68| asymptotically 
attains the form 



{AD-A{0)^D-A{t)) 



SksT {xA- xp) 1 



(70) 



which matches the former expression once to — 
\xA — xd\ and ujo oc 1/ y/\xA — xd\ (see Appendix IE)) . 
We note that, whereas (|68)) does not furnish a pure 
Mittag-Leffler decay, the comparison with the real exper- 
imental data is very good even at short times, as shown 
in [301 ■ Moreover, the expression (pSI) and the corre- 
sponding single-file's [13 provide a compact and elegant 
representation of the D-A distance correlation function, 
reproducing the results obtained formerly starting from 
the same Rouse's protein models [28l[30|. 



D. Correlation function in fluctuating interfaces 

In fluctuating interfaces [8l-ll(j| , often referred as rough 
surfaces, the height h{x,t) obeys to the following gener- 
alized growth equation (generalized elastic line) 



d d"" 

—h{x,t)^-^-^h{x,t)+r]{x,t). (71) 

Accordingly A ix ~ x'j ^ S"^ ix — x' j while z and d are 

left unspecified. The growth of the interface is charach- 
terized by two different growth exponents: the dynamical 
exponent z and the roughness exponent C^. These two 
exponents control the correlations among the neighbor- 
ing heights: the correlation between two different sites x 



and x' grows according to 



t^/^ while the dif- 



ference between the corresponding heights behaves like 



h{x,t)-h(x\t\ 
4J] requires that 



Cr 



A scaling argument 



(72) 



since ea.(|7ip must be invariant under scale transforma- 
tion. However we can obtain both the scaling form and 
the exponents in a rigourous way, by studying the follow- 
ing correlation function 

{[h{x,t)-h{x,0)] h(x',A^h(x',0\ ) = 
{h{x,t)h(x',t'\)~ {h{x,t)h(x',0])- (73) 

^{h(x',A h{x,0)) + {h{x,0) h (x' ,0\). 

Thanks to ([M]) the former function reads (here a = d, 

J — 2z and (3 = {z — d)/z) 



{[h{x,t)-h{x,0)] 
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X — x' 
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(x',0)]) 
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xHf 
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( z-d _\\ 
\ z ' 2z/ 

V 2'2J V~5~'2l/ (-'^^'^'2/ 

X [cos {lo \t — t'\) — cos [ut) — cos {ijjf) + 1] . 

(74) 

a;^/2) and using (IC3|) 



Changing variable I y 



X — X' 



the expression ([71)) can be recasted as 



{[h{x,t)~h{x,0)] h(x',A -h(x',0] ) 
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Jo ^yx 
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(75) 



X cos.[{2yYj^ 
-cos({2yYj^. 



\t-t'\ 



- cos ( {2yY * 
1 
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Eq. (l75t is the generalization of a well-known scaling for- 
mula obtained for the two-point two-time correlation 
function of the Edward- Wilkinson chain [4J, \i^. On 
the other hand the width of a growing surface obeys the 
Family- Vicsek scaling relation [46i] : 



w{L,t) 




dx [h {x, t) ~ {h {x, i))]^ - L'^'-g 



t 

. (^^) 
where L is the maximum size of the system (i"^' is the 

total volume) and g is a scaling function. Comparing 
([75| with ([75]) we find both the scaling expression and 
the correct value of the scaling exponents, i.e Eg. ([7^. 
Moreover we want to stress that, to the knowledge of 
the authors, this is the first time that a scaling function 
appearing in the surface growth is given in an explicit 
form. 

The integrals entering Eqs ([751) and (TM)) can be performed 
with the use of the iJ-functions' properties, and the re- 
sult is expressed in terms of their combination. Here we 
analyze only the simple case of the Edward- Wilkinson 
chain, thus we set z = 2 and d = 1 in ([75| |47| . Thanks 
to (|59l) the correlation function gets the expression 



not only upon time (as the usual FLE do) but also on 
space variable. Indeed, although the FLE reproduces 
the anomalous stochastic motion of the field h{x, t) at a 
given position x disregarding the remaining systems' dy- 
namics, the internal spatial correlations appear through 
the noise term which is correlated in time and in space. 
This is the novelty of our approach: strictly speaking we 
don't loose any information in passing from the Marko- 
vian GEM ^ to the non-Markovian FLE ([21]), on the 
contrary we rather ease the calculations! Indeed, the ap- 
pearance of the Fox function in the two-point two-time 
correlation function, as well as in the noise correlation 
function, makes the computation of any physical observ- 
able relatively easy, if resorting to the few general Fox 
function's properties. Moreover, as in the case of fluctu- 
ating interfaces, the proposed framework allows the ex- 
pression of the statistical quantities in an explicit analyt- 
ical and elegant form involving the iJ-functions. 
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{[h {x, t)~h (a;, 0)] [h (x', t') ~ h (x', 0)]) 



X [cos [Ay'^ 



+00 



dy'- ^ 
-s (By2 



-V2h 



where A 



4#4, 



B 



>in(%/2y) + cos(V22/)] 
- cos (Cj/2) + l] . 

and C = 



(77) 



\x-x'\^' \x-x'Y \x-x'Y' 

Making use of the integrals (|63|) and (|66l) we achieve the 
final scaling form 



{{h (a;, t)~h (a;, 0)] {h (a;', i') - h (a;', 0)]) 
■bJL 

V2 



''''^\x-x'\ 



/(u) = V^e-^-^er/c(y^). 
It is straightforward to verify that 



(78) 



/w 



8-u^/^e ~ M — > 



u — > oo. 



(79) 



VI. 



SUMMARY AND CONCLUDING 
REMARKS 



Appendix A: Hydrodynamic term's Fourier 
transform: case a — d 



If the hydrodynamic term is expressed as 



A(r) 



kt 



(Al) 



its Fourier transform expression reads 



A{q)^{2nr/'\g]'-'/' 



j.d/2 

'^''^d^^Jd/2-1 m r) 



(A2) 

where r stands for \r\ and Jd/2-1 represents the Bessel 
function of order d/2—1. In what follows we will provide 
the exact result for the 1, 2, and 3-dimensional cases. 

i) d — 1. In this case the expression (jA2[) takes the 
simple form 



w N o /""^ , cos(gr) 
Jo a + r 

Solving Ea. (jA3|) we get the following result 



(A3) 



In this article we have shown how the Markovian rep- 
resentation of the system's dynamics furnished by the 
GEM ^ is equivalent to the non-Markovian description 
of the tracer's dynamics given in terms of the FLE (I^TI) . 
Firstly we want stress that the introduced FLE describes 
the time evolution of the randomm field, which depends 



A (q) = 2 sin{qa) ( Si{qa) 1 — cos{qa)Ci{qa) 

"(A4) 
where Si and Ci represent respectively the sine and co- 
sine integrals [3l[ . 
ii) d = 2. We have 
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A (q) = 27T / dr^—^Jo i\q\ r) (A5) 



which gives [42 



Aiq}=2TTKo{\q\a), 



(A6) 



where Kn stands for the modified Bessel function of 0-th 
order [3r|. 

in) d = 3. Eq. (jA2l) takes the following form, 



A(g1 



(27r 



,3/2 



(ir- 



„3/2 



which can be rewritten as 



■Ji/2(k», (A7) 
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sin(|g|r). 



(A8) 



The integral in the previous expression can be split- 
ted into the sum of two contributions, i.e. A{<l} = 



kl 



(/i iq) - h [q)) with 

\h{q) ^aJ-dr^,^£ML 



(A9) 



)(r-x*) ' 



where x — ae"/'^ and x* — ae "/■^. We first study /i (q), 
which can be easily transformed in 



hiq) 



ia\/3 



y 



(AlO) 



where for c.c we denoted the complex conjugated of the 
first integral in the square brackets. After a bit of algebra 
we achieve for the final form of /i (q): 



_ g(a,x) 



gia.,x*) 



g{a,x) ^ f{x)+cos{\q\a)[Si{\q\a)- §]- 
-sm{\q\a)Ci{\q\a). 

(A13) 
We recover immediately the asymptotic expression 

A ((t) ~ ^^^ In (j^^j by expanding Eqs(IS4l),(IMl) and 
the solution of (|A8I) for small |g|. 

Appendix B: Fox function appearance 

The Fox functions are defined as ISSMSa. l4li 



Trra n 
PI 



(ai,Ai) ... {ap,Ap) 

ih,Bi) ... {\,Bq) 



= TT- x{s)y 'ds 



(Bl) 

with l<m<q, 0<n<p. x(s) represents the Mellin 
transform which takes the following form 



X{s) 



UUn+i r (1 - bj - B,s) n^^„+i r (a, + A,s) • 

(B2) 
where Aj and Bj are positive numbers while aj and hj 
are complex. Empty products are interpreted as being 
unity. 

In the expression (|26p the function appearing in the 
denominator can be expressed as a Fox function, see 
ea.p7|. Indeed it is sufficient to notice that the Mellin 
transform of jr-^ is given by 



X{s) 



^,-d. = lifMlii^iM (B3) 
1 + y^ 7 



which matches the definition (IB2p if and only \i m = n = 
p = q = 1, ai = bl = and Ai = Bi = I/7. 



f{x) = cos(|g|a;) [Si{\q]x)^ 



'■] -sm{\q\x)Cii\q\x) 
(All) 
The integral I2 (q) can be manipulated to get 



^2(9) = -^{js^[l-xdy 
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where c.c. this time represents the complex conjugate of 
the whole expression in the curly brackets. According to 
the same procedure used for the former integral /i (q) we 
achieve 



Appendix C: Fox function properties 

In this Section we enumerate the properties of the Fox 
function (jBip that we use throughout our analysis. This 
list is not an exhaustive compendium of the Fox functions 
properties, for which the reader could refer to |34 . l35l . l4l| . 
For convenience in this Section we adopt the following 
short notation 
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[bg.B, 



(CI) 



The useful rules are hereafter listed: 
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The last term in the brackets is transformed with the use 
of den and jCa]), getting finally 



( h{x,Q)-h(x',Q\ ) 
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Expanding the iJ-functions in the former expressions by 
use of (jC7[) we can then get the corrections to 0.5 in the 
brackets. 
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[Up, Ap\ 

KB,] _ 



/ T \ t bi -\-k 
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(C7) 
This expansion is valid whenever J2'i=i Bi ~ X^f^i ^i — 
or ELi 5. - Er=i ^. < and J^j; A, - EL„+i A, + 
sum™^iBi — sumj^^^^Bi > 0. Empty products are in- 
terpreted as being unity. 



Appendix D: Static correlator for fluid membranes 

In this appendix we show how to get a general ex- 
pression for the static correlator (05]). Indeed using the 
property (jC6|) . eq. (|48)) becomes 



( h{x,0)-h{x',0 



1 r+00 dy tt2: 



27r 



{0.25- 



(0,1) (|,^) 

(0,1) (|,i) (0,1) (1,1) 



(Dl) 
Introducing the lower cut-off j/o — x — x' / L and ex- 
tracting the logarithmic term we can recast the previous 
expression as 



Appendix E: Generalized Langevin equation for 

AoA(i) 

In this appendix we derive the generalized Langevin 
equation (GLE) for the single component of the donor- 
acceptor distance /^.u-A^t) in a protein, and then evalu- 
ate its autocorrelation function as arising from the out- 
lined framework. The derivation of the GLE, as well as 
the corresponding correlation function, traces that pro- 
posed in [1^ for single file systems: the main difference 
is that in our analysis we will make use of the Fourier 
transform instead of Laplace in the time domain. 
Subtracting the FLE (PT|) for the donor position h{xD,t) 
from the FLE for the acceptor h{xA, t), the corresponding 
FLE for AD-A{t) is achieved: 



2D'/^AD-A{t) 



CD-A{t) 



(El) 



where we implicitly assumed that xd < xa along the 
protein backbone, and C-D-A(i) = C,{xAit) — C,{xD,i)- 
Eq. (|Eip does not satisfy the generalized FD relation, as it 
is straightforwardly shown by calculating the correlation 
function of the noise 



(Cl?-A(i)Cl?-A(t')) = 

2 {{QixAMixA. t')) - {aXA,mXD,t')}) . 



(E2) 
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The second term in the RHS of the former expression can 
be derived using the general formula ([55)) or by direct 
calculation, once one notices that in this simple case the 
function $(a;,aj) appearing in ^TE\\ is 



^{x,uj) = e 



— „-|x|V-i<^ 



As a matter of fact, in the w space the FD relation reads 



(C(cj)C(w')) = 47rfcsT ^e [-f{uj)] 6{uj + u') 



(E4) 



where 7(w) represents the Fourier transform of the damp- 
ing kernel , which in (|Eip is given by 



(E5) 



Jd-a{^) 



^ -iuj 



and, on the other hand. 



Therefore we sum and subtract on the RHS of (jESp the 
asymptotic expression / ,.^w^ _^ v , obtaining 



(E9) 
(E3) where 7d^-a(^) = 1d-a{'^) - (_,^)(i^,^^) and w^ = 



It is immediate to show that the FD relation (|E4p still 
holds and inverting in time domain we finally get the 
sought form of the GLE for the donor-acceptor distance 
which tends to an FLE with fractional derivative of order 
1/2 in the long time limit fill. [Tgl. l26l[27J 



(ElO) 
We now want to evaluate the donor- acceptor autocorrela- 
tion function: we firstly need the expression of A£i_^(aj) 
from dH), i.e. 



(Cd-a(c^)Ci?-a(^')) = 

8£^^ [1 _ e-xx^-.M (cosxi5-a(w) - sinxD-A(w))] 

(E6) 

with XD^Aii^) = \-\"\-^^'--) _ 

To restore the validity of the FD relation we multiply 
both the terms of (jEip by 



C{^) 



2^ _ ^■s/-iU){Xa-XD) 



2 (l + e-2xD-A(") _ 2e~X"-*(") cosxd-a(w)) 

(E7) 
obtaining the following form of the GLE in the Fourier 
space 

7d-a(w)(-iw)A£,_a('^) = Ci3-a(w) (E8) 

with 7d-a(w) = ^^^(^_^^(.„-.„)-) and Cd-a(w) = 

C(uj)C,D-A(y>)- Now, it is clear that for small uj the 
Fourier transform of the damping kernel gets to the 
asymptotic value 



7d-a(w) 



1 



CD-AH(l-e^^^("»-"°)) 

^D-Aiio) = ^ ^^ '- (Ell) 

\J —ILO 

from whom, thanks to the noise autocorrelation function 

(C^U-Aii^tD-Aiu,')) = 

27rfcBTV2 l-e~^-D-'^'"'(cosxi3-A(c^)-sinxc-A(c^)) (^12) 

|a)|^/2 l+e-2xD-A(")_2e-''-D-A(-) cosxi3-a(w) ' 

we derive the correlation function 

(Az,-aHAz,-a(^')> = ^^3^^ 

X [1 - e-'^"-^^") (cosxd-a(w) - sinxD-A(w))] <5(w + c^') 

(E13) 
Hence inverting in time the previous expression and 

changing variable (y = |w| (a::/i — xd)/'^) we finally 
get 

{/^D-A{t)^D-A{t')) - ^%^(a:A - 2:d)x 



Sr dy 



4y2 t-t' 
'^O'^l 7 7T 



(-iiLj)(a;A -.Tn)' 



l-e-^y (cos \/2y - sin \/2y) 

(E14) 
which exactly matches the expression furnished in (j56p . 
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